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Abstract. Recent developments in off-lattice self-consistent field 
theories for inhomogeneous complex fluids are reviewed. Particu- 
lar emphasis is given to the treatment of intermolecular interac- 
tions and compressibility, to the role of fluctuations, and to the 
discussion of the coarse-graining length which is inherent to the 
theory. Valuable insight can be gained from the comparison of 
self-consistent field calculations with Monte Carlo simulations. Fi- 
nally, some applications of the theory to orientational properties of 
polymers and copolymers at interfaces, and to the phase behavior 
of amphiphiles at surfaces, are presented. 

1 Introduction 

Surfaces and interfaces in supramolccular fluids are a topic of interest in var- 
ious contexts of materials science and soft condensed matter physics. 

On the one hand, the surface structure of materials determines many sur- 
face properties (adhesive and wetting properties, optical properties etc.) and 
thus plays an important role for many applications, such as lubricants, coat- 
ings, thin films and membranes jl| ||. If the topography and chemistry of 
surfaces can be controlled, this can be exploited to taylor surfaces for specific 
purposes[Q, e.g, surfaces which anchor liquid crystals in a well-defined way 
for liquid crystal cells[^, intelligent surfaces^ etc. 

On the other hand, interfacial properties often influence the bulk behav- 
ior of materials in a substantial way, since many commonly used materials 
are inhomogeneous on a microscopic scale. An important class of substances, 
where that is the case, are melts of polymers which contain monomers of 
different type, e.g., polymer blends or copolymer melts. This is because dif- 
ferent organic molecules are often slightly incompatible In homopolymer 
alloys {e.g., composite matrix materials such as rubber toughened plastics), 
the small relative repulsion of the monomers is amplified by the large number 
of monomers in the macromolecules, and completely dominates the entropy 
of mixing, which is proportional to the number of molecules. As a result, 
polymers of different type are usually immiscible at temperatures of practical 
interest. They consist microscopically of numerous finely dispersed droplets 
of one component in the other, and the interfaces between these essentially 
unmixed phases largely govern the material properties @, ||, |^. 

One way to prevent the components from demixing is to chemically bind 
them to each other, thus forming block copolymers. Even though macroscopic 
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phase separation is then inhibited, the chains still have a tendency to rear- 
range themselves so as to allow like monomers to pack next to each other. 
Depending on the temperature, block size etc., this may lead to "microphase 
separation" , i. e, the formation of ordered mesoscopic structures - of lamel- 
lae, ordered micelle arrays, bicontinuous structures jlO| - which can be viewed 
as ordered arrays of interfaces in a sense. The properties of those materials 
strongly depend on the topology and can be quite unusual ||lT|, |2|. Copoly- 
mers are also used as effective compatibilizers in homopolymer blends. Added 
in small amounts, they reduce the interfacial tension between the homopoly- 
mer phases and dynamically prevent the coalescence of droplets (ij, p^ . 
Increasing the copolymer concentration results again in the formation of meso- 
scopically structured microphase-separated phases [T7[| . 

Other prominent examples for inhomogeneous complex fluids are self ag- 
gregating amphiphilic systems, e.g., mixtures of oil, water and soap, or lipid- 
water systems. At high enough amphiphile concentration, they build or- 
dered structures which are very similar to the above mentioned copolymer 
mesophasesJTsI [T9|] . The study of lipids in water is particularly interest- 
ing because of the close connection to biological physics. Indeed, bilayers 
of lipid molecules are a major ingredient of biological membranes and thus 
omnipresent in any living organism. 

Complex fluids are typically made of chemically complicated, large molecules. 
They are usually characterized by a variety of competing length scales - in 
the case of polymer melts the size of the monomers vs. the extension of whole 
molecules, in the case of amphiphilic systems also the correlation within and 
between the self aggregated structures ~ and by the presence of additional, 
conformation and/or orientational, degrees of freedom. A full treatment in 
atomistic detail is in most cases out of reach of today's supercomputers. On 
the other hand, the apparent complexity of the systems actually contributes to 
a simplification of the physics on a coarse-grained level. Due to the large num- 
ber of possible interactions between molecules, microscopic details average out 
to a large extent. A few characteristic attributes of the molecules are often re- 
sponsible for the main features of a substance. For example, the properties of 
polymeric materials are determined mostly by the chain character and the flex- 
ibility of the polymers; amphiphilic systems are characterized by their affinity 
of the amphiphiles to both polar and nonpolar environments, and by their ori- 
entation at interfaces. This motivates the study of idealized simplified models, 
which account only for the main properties of the molecules and absorb the 
microscopic details in a few, effective parameterspO, 21, 2^, ^ |4^. A second 
important point is that dense macromolecular systems are often unusually 
well described by mean field approximations. Since large molecules interact 
with many others, the effective interaction range in the limit of high molecular 
weight is very large, and according to a simple Ginzburg type argument, the 
critical region in which concentration fluctuations become important is very 
small as a result [p3| . 
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Mean field theories are thus widely used to describe inhomogeneous macro- 
molecular systems. Many simple and largely successful theories are based on 
Ginzburg-Landau type functionals of the concentration or other extensive 
quantities (fields), which characterize a system with respect to the problem 
of interest [^. The interactions between molecules are effectively incorpo- 
rated into square gradient terms and sometimes contributions of higher order 
derivatives of the fields The parameters of the model are either derived 
from a microscopic theory (as in the de Gennes-Flory-Huggins theory of in- 
homogeneous polymer mixtures or from comparison with experiment. 
Such a treatment is adequate as long as the interactions are short ranged, 
the inhomogeneities are only weak and characterized by length scales of the 
order of the molecular size. It becomes questionable in systems which are 
inhomogeneous on smaller length scales - which is true for most examples 
cited above. A substantial amount of effort has therefore been put into de- 
veloping more refined mean field theories, which explicitly account for the 
chain conformational structure (lattice theories such as the Flory-Huggins 
theory [psf pO| , p6| , |27t , density functional theories p^, [2£ , 30, Bl], self-consistent 
field theories [p2[ |3^ , |3^ , p6|), or even for the local structure in polymer 
melts (e.g., the lattice cluster theory by Freed and coworkers |37t , or the P- 
RISM theory by Schweizer and Curro[Q). The more microscopic factors a 
theory incorporates, the more detailed is the information that it can provide 
on local structure properties. Unfortunately, the treatment also gets more 
and more involved. Theories with a higher level of coarse-graining have the 
advantage of being somewhat more transparent and easier to handle. More- 
over, the origins of given physical phenomena can sometimes be tracked down 
more easily, if they can be discussed in terms of a highly idealized model. 
The optimal choice of a theory of course depends on the specific questions 
one wishes to address. 

In this paper, we shall discuss one of the most powerful mean field tools 
in the study of macromolecular systems, the self-consistent field approach. 
Even though, strictly speaking, every mean field approximation involves self- 
consistent fields, the name of self-consistent field (SCF) theory is in the poly- 
mer community by tradition reserved to a certain type of approach, which 
goes back some time ago to Edwards |^ and Helfand and Tagamifs^. In a 
nutshell, it can be characterized as follows: Polymers are described as random 
walks in a positional dependent chemical potential, which depends in turn in 
some self-consistent way on the distribution of monomers. The basic idea is 
thus extremely simple, which explains the enormous appeal that it has had for 
polymer scientists ever since it was first formulated. A particularly popular 
version of the theory has been developed by Scheutjens and Fleer for lattice 
models [p6| and applied to various problems related to polymers at interfaces 
or surfaces |^ . Lattice models allow for an efficient study of systems with a 
planar geometry, as long as one is not crucially interested in chain stiffness and 
local orientation effects. The treatment of situations which involve strongly 
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curved interfaces and local or global orientational order is more difficult. One 
way out of this dilemma is to employ more sophisticated lattice models. For 
example, the bond fluctuation model of Carmesin and Kremerp3[, a popular 
lattice model in Monte Carlo simulations! 4L has recently been used for 
self-consistent field studies of lipid bilayers and polymer alloy systems |Q. On 
the other hand, the above mentioned difficulties are automatically eliminated 
if one works in continuous space. We shall restrict our discussion mostly to 
such off-lattice models in the following. The basic concepts and recent me- 
thodical advances shall be reviewed in the next section. In section 3, some 
applications shall be presented which illustrate the use of the model for the 
study of polymer interfaces and amphiphilic systems. We summarize and 
conclude in section 4. 



2 Theory 

2.1 Basic concepts 

We consider a mixture of rij molecules of type j = a,b,c - ■ ■ in a volume 
which are built from Nj monomers of species a = A, B,C ■ ■ ■. Here ro- 
man indices distinguish between molecule types, and greek indices denote 
monomer species. For simplicity, we shall also assume that the monomers 
in a molecule build a linear sequence ~ the generalization to different archi- 
tectures is straightforward . The molecules can then be represented by 
continuous space curves R{s), with s varying between and Nj. Polymer 
types J differ from each other by their chain length Nj^ and by the distri- 
butions of monomers a on the chains. The latter are conveniently described 
by functions ^a.j{s), which are 1 for portions of the polymer occupied by the 
a-monomers, and otherwise. Hence one has J2a^a,jis) = 1 for all j and s. 
In A-homopolymers, for example, 7a(s) = Sa^o for all s € [0,Nj], and in the 
case of symmetric A : B diblock copolymers, 7q(s) = 5A,a for s e [0, iVj/2] 
and 7q(s) = 5B,a for s & [Nj/2, 1]. 

For a given configuration {Rij{-)}, one can then define a local monomer 
density operator 

Pair) =J2^^ P'o^^^ (1) 

3 

with pi{f) = — 22 ^(^^ 7a,i(s)- 

The interaction between monomers can be described by a coarse-grained 
monomer free energy functional V{pa}- In addition to the Boltzmann factor 
associated with this interaction energy, chain conformations are distributed 
according to an intrinsic statistical weight which accounts for the 

internal energy of the chain, and for the part of the configurational entropy 
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which steins from length scales smaller than the coarse-graining length. In 
practice, polymers are most commonly modelled as gaussian chains 

P,W)}=AAexp[-5]^ pds\^\\^,is)], (2) 

with the statistical segment length of a-monomers ba, and the normaliza- 
tion factor A/" chosen such that / 'D{R{-)} Vj{R{-)} — V. Here and in the 

1/3 

following, all lengths are given in units of some microscopic length Vq and 
hence dimensionless. The conformational free energy of the chains is thus 
approximated by a gaussian stretching energy with spring constant 3/26q,^. 
Shorter or stiffer polymers are sometimes represented as "wormlike" chains of 
fixed contour length. Chain portions made of a monomers are then charac- 
terized by the monomer length Cq, and the dimensionless stiffness r]a , and the 
probability distribution functional Vj are given byji^ 

s a ^ 

where U — {dR/ds)/aa is a dimensional tangent vector constrained to unity 
by the delta function. This choice of a statistical weight has the advantage 
that monomer orientations are well-defined and orientation dependent inter- 
actions between monomers can be introduced in a straightforward way. In the 
limit of small rj, i.e., very flexible chains, the wormlike or semiflexible chain 
statistics reduces to gaussian chain statistics[^. Of course, it is also possible 
to consider chains with mixed gaussian and wormlike parts. 

Alternatively, many polymer models describe molecules by discrete walks 
{Rs} instead of continuous curves {R{s)}. The treatment is essentially the 
same, save that integrals J^^ ds are replaced by sums J2s ■ '^^ freely 

jointed chains with fixed bond length Z, for example, the statistical weight is 

V,{R{-)}^X II s{\R,+,^Rs\-l). (4) 

s=l 

The introduction of bending potentials acting on the angle between subse- 
quent bonds, monomer dependent bond lengths laf3 (for bonds connecting 
monomers a and f3), or variable bond lengths with bond length potentials, 
etc., is straightforward. 

In all the examples quoted so far, the chain statistics is determined by 
local relations between neighbor monomers along the chain, and is thus that 
of random walks. As we shall see in section 2.2, solving such a model in 
self-consistent field approximation amounts to solving a diffusion equation 
self-consistently, which can usually be done at relatively moderate expense. 
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This makes this kind of approach particularly attractive. However, the self- 
consistent field method is not restricted to locally defined weight functions. 
Szleifer and coworkers |5^ have recently devised a procedure which allows 
to perform self-consistent field calculations for chains with arbitrary chain 
statistics, e.g., self avoiding walk statistics etc. (see section 2.2). 

In the following, we shall use the notation JVj{R{-)} = JV{R{-)}'Pj{R{-)} 
for the weighted sum over all chain conformations, which is a path integral in 
the case of continuous curves, and a summation over discrete walks otherwise. 
The partition function of the system in the canonical ensemble is then given 

by 

71 j 

exp[-/3V{p4], (5) 
where (3 = l/ksT is the Boltzmann factor, and the monomer interaction func- 



tional V{pa} has been introduced earher. In order to proceed, it is useful |34 
to insert functional integrals over delta functions 

1 = / V{p^} 5{p^ -p^)= I V{p^} [ P{c.J e^^^-°« (P=«-P.(.l)^ (g) 



which allow to replace the density operators Paif) in V{pa} by density func- 
tions Pa (r) . The subscript indicates that the limits of integration of the 
auxiliary fields LUa{r) are —ioo to ioo. Eqn (|^) can then be rewritten exactly 
as 

2c = { n / ['^M} exp\~ f3J^c{Pc.,u;a}] (7) 
with the canonical free energy functional 

f3J^c{Pa,^a} =^ PV{Pa} - J df IVair) Pa{r) - Hj ln{Qj {uja} / rij) . (8) 

a j 

The functional Qj{u!a} is the partition fimction of a single chain moving in 
the (imaginary) external fields LUa{r), 

QjiiVa} = J 'Dj{R{-)} exp 

'rfsco„(i?(s))7„j(s)]. (9) 

The main step of the self-consistent field approach consists in performing 
a saddle point integration of the integral (Q) with respect to oja (r) : the path 
integral is replaced by the value of the integrand at the corresponding saddle 
function ujai'r') in the complex plane. Extremization of the exponent in (|^) 
yields the equations 

Sln{Qj{uJa}) 1^ ,-,„^ 
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Brackets {■)c denote statistical averages taken in a canonical ensemble of Uj 
noninteracting chains of type j subject to the external fields Wq (r) . Note that 
the saddle functions LUai^r) which solve these equations are now real. The 
approximation thus amounts to replacing the exact constraint Pa{r) = Pa{^ 
in eqn. by the more relaxed requirement = {pa{r))c- The functions 

LOa (r) act as Lagrange parameters which enforce this condition. The resulting 
partition function is 

2c = { n / '^iP^} } exp [ - /3 .Fc{pa, uj^Ipo.}} ] , (11) 



where the fields uja{pa} depend self-consistently on pa according to eqn. (|10|). 

When looking at mixtures, it is sometimes more convenient to work in 
the grand canonical ensemble p5|, p2[ . The starting point is then the grand 
canonical partition function 

^Gc-n E l^rr n /^.{^a-)}} eM-mpa}], (12) 

where fj,j are the chemical potentials associated with molecules of type j. One 
can now proceed the same way as above, and obtains the analogue to (|^) with 
the grand canonical free energy functional 

a j 

The extremization with respect to LUa{f) yields the condition 

p^if)^-Y^e'^^^^ = {Mr))ac ■ (14) 

Hence the grand canonical partition function in self-consistent field approxi- 
mation reads 

-ZgC = I /^?{/0q} I exp - P TGc{Pa,^^a{Pa}} ■ (15) 

with uja{pa} defined via eqn. (p^). 

At this point, concentration fluctuations are still included in the partition 
function. Usually, self-consistent field approximations take one more step and 
perform a second saddle point integration of ( [Til ) or (15) with respect to Pa{r)- 
The free energy is then approximated by 

F^~kBT\iiZKF^^^ = min.F{p„,Wa{Pa}}, (16) 

Po(r) 
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with the minimization equations 



^Jr) = J — T-T (17) 

in both the canonical and the grand canonical ensemble. This is a very in- 
tuitive, typical mean field result: The auxiliary fields cjQ(r) which drive the 
monomer densities pa (r) are identified with the local excess free energy needed 
to add one monomer a to the melt at the position r, less the contribution of 
the translational entropy of polymers. Since the fields LUa are functions of pa 
by means of eqn. ( |lO| ) or (E4), eqn. (17) effectively determines the mean field 



concentration profiles Pa'~^\f). 

Hence two approximations enter the self-consistent field theory: The first 
leading to eqn. ( [Til ) or eqn. (p^), respectively, and the second leading to eqn. 
(©). We have taken the care to separate these two steps for the following 
reasons. First, Shi, Noolandi and Desai have recently suggested a way how 
to systematically improve on the second step and take Gaussian or even higher 
order concentration fluctuations into account (see section 2.4). This allows 
for a stability analysis of the mean field solution|Q, and for the construction 
of kinetic paths connecting different mean field solutions ||55[|. Second, the two 
steps differ qualitatively from each other. Technically, the saddle functional 
approximates integrals of lo^ on the imaginary axis, and real integrals of Pa. 
A stable saddle point is thus a minimum with respect to Pa and a maximum 
with respect to uja- The first step motivates the introduction of a coarse- 
grained density functional T{paTUJa{pa}}i which is then minimized in step 
two as in usual mean field approaches. 

The introduction of orientation dependent monomer interactions into the 
formalism is straightforward. The monomer density operator (0) is simply 
replaced by 



1 fN, 

pi{r,u) = — Y^ ds,5(f-4(s))<5(u-C/,^.(s))7„,,(s), (18) 
where U is the orientation of monomer ij. The fields uja become orientation 



dependent and meet (c/. <\l7\)) uja{r,u) = S (3V / S pa{r, u) . 

Finally in this section, we sketch Helfand's original formulation of the 
self-consistent field theory Jsst, which is particularly suited for the study of 
polymer interfaces and surfaces. Consider an interface between two entirely 
segregated homopolymer phases j = A and B, which occupy each a partial 
volume Vj = rijNj/p*. The two pure systems are characterized by their bulk 
densities p* and the associated fields lu* — dv*/dp*, where v* is the bulk 
density of the interaction energy V. Let further the probability distribution 
'Pj{R{-y\ be normalized such that J I?j{/i(-)} = Vj. The free energy density 
of the pure systems is then given by /?/* = (iv* + p* /Njln{p* /Nj). When 
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looking at the interface between the two phases, it is now convenient to work 
with shifted monomer fields 



^3{r)=^j{r)-uj], (19) 

which vanish in the pure system j . They can be interpreted as the excess work 
needed to bring a monomer of type j from the pure system into the interfacial 
region. The underlying shifted interaction potential {oJj{r) = 5V / 5pj{r)) has 
the form 

/3V = /3V - ^ uj* j drpj (f) - const., (20) 

3 

where "const." = Vj { [3v* — } subtracts the bulk contribution of the 
pure phases. The shifted fields ujj define a renormalized single chain partition 
function 

Q, = j V,{R[-)} exp [ - j\s ti,(i?(s))] = Qj exp[-a;;iV,], (21) 

Note that the actual value of Qj is Qj — rijNj/p* Hence the equation for the 
density profiles can be written in the simple form 

P3irl = -^^^ (22) 
and the excess free energy of the interface is 

The formalism can also be applied if the polymers A and B are not fully 
segregated, i.e., the bulk phases contain both A and B. In that case, the 



interfacial free energy (23) has to be corrected for additional bulk contribu- 
tions. Helfand et al have approximated the excess interaction potential V by 
the simple form^, ^ 

PV = / dr + ^(^ + ^ - 1)^ }, (24) 

J \JPaPb ^i^^bT pg 

where the Flory-Huggins parameter x describes for the incompatibility of 
monomers A and B, and k is the compressibility of the melt. In the limit of 
infinitely long polymers and zero compressibility, the theory could be solved 
analytically by Helfand and Tagami (SSL: strong segregation limit ]3^ ) . For a 
symmetric mixture (6a = — b, p\ ^ p*g ~ p), they obtained the interfacial 
tension Pctssl — ■\/x/6 pb and the interfacial width wssl = b/^/Gx. 
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2.2 Implementation 

Eqns. @ or dlil ) and (|l7| ) define a complete cycle of self-consistent 
equations, which can be solved numerically by suitable iteration methods. The 
main task in each iteration step consists in evaluating the partition function 
Qj{uja} of a single chain j in the external fields uJa{r) (^) and the derivatives 
of Qj{uja}- If the polymers have random walk statistics VjiRi )} as in (||), 
(H) or it is useful to define the propagators 

G,{r,t;/,t') = ypj{i?(.)}exp[-^^ dsc^<,(i?(s))7„,,0 

X 5{f- R{t)) Sir' ~ R{t')), (25) 
or, in case monomer orientations are important, 

X S{f- R{t)) Sir' - Rit')) Siu - C/(t)) <5(u' - C/(t')). (26) 

Since they carry no memory along the chain, these propagators satisfy mod- 
ified diffusion equations with the initial condition Gjir^t^r'^t) = Sir — 7") or 
Gj(f, M, t] f" , u' , t') = Sir — 7") Siu — u'). Specifically, the diffusion equation 
for Gaussian chains (||) reads |Q 

+E^"^^(*')[- f +^"(^")]] G,(r,t;r',t') = 0; (27) 
for wormlikc chains 

^ + E K^'Vr' - ^V|, +cc>4r^)]] G, (f, u, t- r^, u' , t') = 0, (28) 

a 

and for freely jointed chains with bond length I (||), 

Gj(f,<;r',t' + 1) = cxp[-tJc«(r')] f du GAf,t;r' - I u,t'), (29) 

where |m| = 1 as before, and J du denotes integration over the full solid an- 
gle. Once the propagators have been determined, one easily calculates the 
partition function Qj — Jdrdf" Gj(r, 0; r^, iVj) and its derivatives 

SQ- f f^^ 

dfdr' I dsjajis) Gjif,0;fo,s) Gjifo,s;f',Nj) 



Suiairo) 

^ - 2 drdf* dsjajis) j ds' -ip.jis') (30) 



SuJairo)5u:piri) 



Gjir,0-ro,s') Gj(fo, s'; fi, s) Gjin, s-y , Nj) 
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etc. In fact, only the first derivatives enter the self-consistent field equations. 
It is thus often less time consuming to calculate directly the "end segment 
distributions" qj{r,t) = Jdf* G j {r' , 0; f, t) a.nd qj{f,t) = Jdr'Gj{f,t;f',Nj), 
which solve the same diffusion equation as G with boundary conditions qj (f, 0) = 1 
and qj (r , Nj) = 1 . The end segment distribution functions can be used to 
calculate density profiles, chain end distributions etc. The full propagators 
Gj{f,t;r',t') are needed when looking at the effect of fluctuations (see sec- 
tion 2.4), or at chain correlation functions such as distributions of end-to-end 
vectors. 

Diffusion equations can usually be solved without too much computational 
effort. The computer time needed for one iteration step only scales linearly 
with the chain length N , hence one can handle relatively long chains. Physi- 
cally, treating the chains as random walks is well motivated in dense melts of 
long polymers, since the excluded volume interactions between monomers of 
the same chain are screened by the presence of other chains 

However, this is not true for dilute polymer chains, and on length scales 
smaller than the screening length||2l|. The random walk approximation thus 
becomes questionable when looking at polymers in good solvent conditions, 
or at relatively short macromolecules. In that case, it is more appropriate to 
choose a weight distribution Vj which accounts also for long range correlations 
within a chain. One then has to consider whole chains in the external field of 
the other chains (single chain mean field). 

Szleifer and Coworkers have recently suggested an enumeration procedure 
for the evaluation of chain partition functions |5^, ^ . They approximate the 
path integral Jl?j{i?(-)} by the sum over a representative sample of chain 
configurations Rij{-), which are distributed according to the weight function 
7'j{i?(-)}. The sample can be generated by Monte Carlo simulations, taken 
from experiments etc. This approach is conceptually extremely simple and 
can be applied to arbitrary weight distributions Vj{R{-)}. For example, dilute 
polymers in good solvent are well described by self avoiding walks, and one can 
draw the sample from single chain simulations pO|| . Semidilute polymers have 
self avoiding walk statistics on short length scales below the screening length, 
and random walk statistics on larger length scales. The best description on 
all length scales is obtained when the sample is generated in Monte Carlo 
simulations of a melt |46| . Unfortunately, the computational effort required 
to reliably solve the self-consistent equations grows rapidly with the chain 
length, since the sizes of the sample have to be made very large. In dense 
melt of long polymers, random walk models for chains usually give reasonably 
good results at much less expense. 

Nevertheless, the method is very versatile and has much potential for fur- 
ther refinement. For example, reweighting schemes can be devised in analogy 
to the reweighting techniques in Monte Carlo simulations |Q, which are use- 
ful if the actual distribution of chain conformations differs very much from 
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the distribution 'Pj{R{-)}. One then draws the sample of chain conformations 
according to a modified distribution V' j{R{-)}, and corrects for this in the 
sum over conformations. 

jv,{R{-)} — ••• 

Weinhold et al have introduced a couphng between Vj and the number of 
interchain contacts, and could qualitatively reproduce the shrinking of chains 
in an athermal melt with increasing monomer density ]5^. Within a similar 
approach, Miiller calculates the statistical mechanics of clusters of polymers 
embedded in a self-consistent field, in order to study the conformations of a 
single A chain in an unfavorable B environment [^Sj . His results compare well 
with Monte Carlo simulations. 

With these remarks, we end the discussion of the Szleifer method. The 
reader is referred to ||5l| for a detailed review of the approach and its appli- 
cations. We close this section with a few technical comments on the solution 
of the mean field equations. 

The careful choice of the suitable basis is crucial for the success of a 
self-consistent field calculation. Matsen et al have implemented with great 
success a method to study periodically ordered mesophases, where spatially- 
dependent functions are expanded in orthonormal basis functions with the 
symmetry of the phase being considered p9|, |60t |6l| . If monomer interactions 
come into play, it is usually beneficial to expand functions of the orientation 
u in spherical harmonics ||6^. 

The mean field equations are solved iteratively, e.g.^ using the fields uja{f) 
as iteration variables. The most popular iteration procedure is the Newton 
Raphson method One usually obtains a solution of reasonably high ac- 
curacy within less than 10 iteration steps. However, each step involves the 
determination of a full Hesse matrix of derivatives, i.e., the number of func- 
tion evaluations grows quadratically with the number of degrees of freedom. 
Relaxation procedures take more iteration steps, but may nevertheless require 
fewer function evaluations in the end. The present author has had good expe- 
rience with a variant of a method originally suggested by NgQ, Another 
efficient procedure is the multidimensional secant method of Broyden||66|, |67| . 

2.3 Monomer interactions and compressibility 

We now turn to the discussion of the functional V{pa}, which defines 
together with 'Pj{R{-)} the actual model. It describes the free energy of a 
system of interacting monomers without translational entropy. Naively, one 
could be tempted to identify it with the sum of direct interactions Wapif) 
between monomers, ^ Jdrdf" J2a/3 ^a/3{r — f')pa{r)pp{f'). However, this ap- 
proximation is poor in dense fiuids, where indirect interactions and multiplet 
correlations are important. In most applications of the self-consistent field 
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theory, the quantities of interest are composition inhomogeneities, whereas 
the total density varies comparatively little throughout the system. Further- 
more, the relative repulsion of monomers of different type is usually vanish- 
ingly small compared to the total free energy in the melt. It is thus propitious 
to separate V into an "equation of state" contribution Vo, which accounts for 
density fluctuations and compressibility effects, and an interaction part Vi„t„„ 
which incorporates the incompatibility of unlike monomers, 

V{pc.} = Vo{p„} + V.„,„{p,}. (31) 

We shall discuss both parts in turn, and begin with Vq. If the melt is 
nearly incompressible and the density fluctuations are small, one may use 
Helfand's quadratic approximation |3|, ^ (see eqn. (|2^)) 

Vo{pa} = ^ /"dr (l-^p„z;„)2, (32) 

where Wq, is the specific volume of an a monomer. 

More rigorously, total incomprcssibility is often postulated. 

'^^Pai'r) = 1 everywhere. (33) 

This is most commonly achieved by introducing an additional Lagrange pa- 
rameter field £,{r), which couples to the incomprcssibility constraint (|33|)[p4| 
^ p9| . The free energy functional (||) or ( ll3| ) is then replaced by 

= /3V.„,„{pJ - /dr a;„(r)p,(f) - /df ^(f) [E„ P^l^) «a - l] 
'^jnjhi{Qj{uja\ /rij) (canonical ensemble) 
X]j e'^^J Qj{(jjQ,}. (grand canonical ensemble), 

and the minimization equation (|l7| ) turns into 

Wa(r) = '^/"'^rlr" - ^{r) Va- (35) 
dpa[r) 



The incomprcssibility constraint (^3|) determines the field ^(r) unambiguously, 
except for a constant ^o- It is often chosen such that / dr£,{f) = 0. Note that 
the choice of determines the offset of the chemical potentials pj in the 
grand canonical ensemble. 

Within this framework of an incompressible theory, total density fluctu- 
ations is often reintroduced by means of a hypothetical noninteracting "sol- 
vent" S, which fills the space between the molecules p^, |5l| , |68| , |6^ . The 
incompressibility condition (^ ) then reads J2a P'^i''^ '"<^ + Ps{r)vs = 1, and 
the excess potential of the solvent is simply given by lOs (r) = ^(r) Ws • Ac- 



cording to eqn. (10) or (Qj), the density of the solvent is thus given by 
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Ps{r) = z exp[— ^(r) Vs], with z — Us/ Jdrex.p[—£^{f)vs] in the canonical case 
and z = explPus] in the grand canonical case. This defines ^(r) as a func- 
tion of the monomer densities Pair). Choosing i^o such that ln{zvs) = 0, one 
obtains the fields 

^c{r) = ^.f^rT ~ - Ea ^«)- (36) 

The approach amounts to using an equation of state potential of the form 

/3Vo{pJ = - /dr ( 1 - r;(r) ) In ( 1 - 77(r) ) + - / drr,{r}, (37) 

Vs J Vs J 

where the packing fraction ri[f) ~ Pa{'r)va is the local volume fraction oc- 
cupied by monomers. The last term ensures that the free energy per monomer 
vanishes in the limit 77 — > 0. Eqn. ( |3^ ) is a local version of the familiar 
compressible Flory-Huggins or Sanchez- Lacombe theory [f70|. It provides an 
intuitive and straightforward treatment of compressibility effects in polymer 
blends. Unfortunately, the prediction for the equation of state of the melt 
can be rather poorjT^, due to the total neglect of details of the monomer 
structure and local monomer correlations. 

An alternative approach suggested in Ref. |Q conceives Vq as the den- 
sity functional of a suitable reference system of identical monomers with no 
translational entropy. In local density approximation, Vq can be derived from 
the equation of state 11 (p) (11 is the pressure), 

/3Vo = [drp{r) f[p{r)] with J{p) = ("dx H^. (38) 







The function f{p) is the local free energy per monomer. When studying 
interfaces or surfaces, Vq is conveniently replaced according to ( pO| ) by 

/3Vo = jdf[ p{r) ( f[p{r)] - lo*) - const.] , (39) 

where uj* = f{p*) + p*{df/dp)p-' is the excess free energy per monomer at 
the bulk density /?*, and const. = p*^(df /dp)pt . The theory requires as input 
the knowledge of the equation of state of the reference system, i.e., a melt of 
infinitely long polymers built from the reference monomers. For example, one 
can take the Flory-Huggins form|^ 

Hip) Ip=-\- l/v ln(l - 77), (40) 

with the packing fraction 77 — p/p*, and recover the potential ( ^7| ) with 
Vs ~ Based on the Carnahan-Starling equation of state for single 

hard spheres, Dickman and Hall have derived an equation of state for hard 
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chains[^, which compares well with computer simulations of hard chains, 
and even of the bond fluctuation model ifzil], 



n{p) /p - Co 



1 + + r/ — ?7 



(41) 



and leads to the free energy 

/(p) = Cor;(4-3r?)/(l-r;)2. (42) 

The constant Co depends on details of the intramolecular chain structure, 
and rj = a^pn/G is the actual volume occupied by monomers of diameter a. 
Other forms for the equation of state are available as well[Q. In some cases, 
one may wish to include more information on the local liquid structure in 
the functional Vo, and in the corresponding monomer excess free energy w(r). 



Nath et al 75 1 suggest a modified version of the self-consistent field theory. 



where the excess free energy is given by (cf. ( |19| ) 

Cj{r)^- fdr' c{r-r')[p{f)- p*]+ij*, (43) 



with the direct correlation functions c{r) taken from P-RISM theoryPq]. 

Next we discuss the second contribution Vi„ter to the free energy functional 
V{pa}, which describes the monomer specific part of the interactions. Since 
it is small compared to the total free energy of the system, we may assume 
a perturbative treatment. We split the pair interactions Waf^ir) between 
monomers a and /3 into two parts Wafs (r) = Wq (r) + W^p (f) . The potential 
Wo{r) describes the interactions in the reference system Vq. Then we define 

/3 C/a/3(r) = 1 - exp[-/? <^(f)]. (44) 

If W^f^{'r) is integrable and small, Uoip{r) reduces to VF^^(r). More generally, 
the expression (|44| ) also allows to deal with effects of non-additive packing, 
size disparities of monomers and the like (see also Ref. ^ij) In perturbation 
theory the free energy contribution Vi„tc, is given by 

V,„.„{p4 = i / drd^ ^) /'S(^^' (45) 

(2) ... ^ 

Here p)^^ is the pair density of type a monomers at r and type /3 monomers 
at r*, which are not direct neighbors along one polymer chain. (The interac- 
tions between the latter contribute to the conformational weight functional 
'Pj{R{-)}). Furthermore, the pair distribution function is approximated by 

^) = P-(r) pp{^) l{r- r'), (46) 
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i.e., the monomer pair correlation fmiction 7(r) is taken to be independent of 
the identity of the interacting monomers, and of their densities ||77|]. 

If the interactions are short range, the profiles ppir') can be expanded 
around r.[03| With the definitions 



Xa/3 = (3p* jdrUo.p{^-i{r) (47) 
<0 = Pp* !drU^p{r)-i{r)nr, (48) 



etc., one obtains 



The indices i run over cartesian coordinates x, y, z, and the bulk density p* 
was introduced in order to make Xa[3 dimensionless. The resulting excess 
potentials are given by 

ij •' 

In systems with only two types of monomers A and B, the parameters Xq/3 
are usually combined to one single "Flory-Huggins parameter" 

X = XAB -'^{xaa + Xbb)- (51) 

Similarly, we define cr*-' — ct^^ — ^i'^AA + "^bb)- "^^"^ fields o^a.b can then be 
rewritten as 

= ^ + :Jr*PB + 5irE--'S!: + '-W0} (52) 



(5p p* 2/9* -^--^ dridrj 

ij 



where fe{p(f)} = ^(x^A - Xbb) P + ^ ^(g^^ - ergg) J (53) 

depends only on the total density profile p{r) = pA + ps, and contributions 
which are identical for both components have been dropped except for the 
leading term 5(3Vt^/5p. Similarly, the free energy functional Vi„tor reads 

/3V.„.„ = ^ ( dr[x PAPB-\Y.<r^''^'^+h{p} {pA~PB)+v{p}]. (54) 
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Here again, v{p} can be neglected compared to /3Vo. 

Usually, only the "local" contributions (cx Xaff) are taken into account. 
However, the higher order nonlocal terms become important in situations 
where the local density profiles vary strongly. We shall illustrate this with 
the example of surface segregation in a completely symmetric binary {A, B) 
blend of two incompatible, but otherwise identical homopolymers. Even if 
the surface is chosen neutral, i.e., the surface tensions of A and B are equal, 
the minority component segregates to the surface, because it has less unfa- 
vorable contacts with polymers of the majority component there. This effect 
is often called a "missing neighbor effect" . It increases with the range of the 
interactions, i.e., it should increase with the value of a. 

In Ref. self-consistent field calculations are presented for such a sys- 
tem. The parameters of the theory were adjusted to the bond fluctuation 
model, a lattice model for polymers which has been widely used in Monte 
Carlo simulations and has extremely well known bulk properties. It will be 
introduced in some more detail in Section 3.2. The chain conformations in 
the bulk[f78|, the equation of state[^ the thermodynamics of mixing 
have been investigated in detail. Therefore, the model parameters for the 
self-consistent theory, e.g., within the Gaussian chain model, are all known. 
Fig. |l| shows two segregation profiles at bulk two-phase coexistence, one with 
(T = (dashed line) and one with cr adapted to the bond fluctuation model 
(solid line). (Note that the field h{p} (|53|), which couples directly to the com- 
position fluctuations, vanishes in a perfectly symmetric mixture.) In the limit 
of tr = or pure contact interactions, surface segregation is almost entirely 
suppressed. Using the adjusted value of cr, one finds that the volume fraction 
'^A — PaI p of the minority phase increases by almost a factor of two at the 
surface. 

The self-consistent field results can be compared with Monte Carlo simu- 
lations of Rouault et aZj8|]. Fig. || shows the predicted difference between the 
volume fractions of A at the surface and in the bulk A$yi, and corresponding 
Monte Carlo data. The simulations were performed in a slab geometry with 
a slab thickness of about three times the gyration radius of the chains. In 
such a thin film, the composition variable \pA — Pb\ at coexistence is reduced 
compared to the bulk system. The difference between the volume fraction <I>^ 
at the surface of the thin film, and in the bulk of the infinite system, thus 
gives an upper bound for the actual value of A$^. The difference between 
4>yi at the surface and in the center of the thin film gives a lower bound. 
Fig. H illustrates that the theoretical prediction lies nicely within these two 
bounds. We emphasize that the good agreement was reached without any 
adjustable parameters. However, the range of the monomer interactions has 
to be accounted for correctly, i.e., the nonlocal term a may not be neglected 
in this situation with strong density variations. 

A second caveat applies to the approximation (^6|). In general, the pair 
correlation function depends on the identity of the monomers, and on the 
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local environment. For example, the number of monomer contacts divided 
by the local density typically decreases in the vicinity of an interface [Q, 
and is different for homopolymers and copolymers [^Sj. This is a higher order 
effect in our perturbative approach, which is usually not taken into account 
in self-consistent field theories. Wayne et al have recently found experimental 
evidence that it may engender a difference of the effective x parameter in 
homopolymcr and block copolymer blends [Q. 



2.4 Fluctuations 

In the limit of infinite chain length (while keeping the densities constant), the 
saddle point integrations leading to the free energy expression ( p^ become 
exact. At any finite chain length, however, the free energy is affected by con- 
centration fluctuations. On the one hand, the concentrations always fluctuate 
locally about their local mean value. The length scale for these fluctuations 
is given by the bulk correlation length, which is usually roughly the gyration 
radius and diverges only very close to a critical point. In systems with in- 
terfaces, a second qualitatively different type of fluctuations emerges: Since 
the interface breaks a continuous symmetry, the translational symmetry, long 
wavelength goldstone excitations come into existence which cost virtually no 
energy in the zero wavelength limit. The corresponding length scale thus di- 
verges at all temperatures. These fluctuations are sometimes referred to as 
interface "waviness" , or capillary wave fluctuations. 

Local concentration fluctuations can be assessed systematically by an ex- 
pansion of local concentration perturbations, 5pa{r) = Pq(^) — pf^^C^), about 
the SCF concentration proflles pf^^ {r) . A way how to perform such a pertur- 
bative extension of the self-consistent field theory has recently been suggested 
for incompressible mixtures in the canonical ensemble by Shi, Noolandi and 
Desai||5^. In the following, we shall briefly sketch their approach, and general- 
ize it for the case of compressible mixtures and arbitrary statistical ensemble. 
Capillary wave fluctuations shall be discussed subsequently. 

The starting point is the partition function ([l^) or (^sj). In the perturba- 
tive treatment, the free energy functional T{pa,0Ja{pa}} is expanded about 
its self-consistent field minimum F^^^ . Following Shi et al, this is done most 
conveniently by expanding !Fc{pa,^a} in both 5pa and the field deviations 
5uja = ^a{Pa} — ^a{Pa'^^}, keeping in mind that the latter are functions of 
the concentration perturbations Spa- It is useful to define the single-chain 
fc-point distribution fimction 

Sai-.-CfcC^l • • -n:) = 7^ T ,^ f ^'l 7^, (55) 

which is the joint fc-point density (p^j(n.) • • -pli^i^k)) for monomers of the 
same chain of type j. The corresponding cumulant correlation functions are 
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given by 

^ 1^ (Kf^\ 
Cii-a^l'^l ■ • -^fc = T Fm X Fm- 

For example, the lowest cumulant correlation function is simply the den- 
sity, c^(r) = gii{r) = {pii)j the pair correlation function is c'^p{'r,r') = 
^ etc. Note that these are single chain correla- 

tion functions of noninteracting chains j subject to the external fields ^^(r). 
For convenience, we also introduce the notation 

'^ai---ak [ri-'-rk)^ ,^ / ^ T=rT (57) 



and define the functions Ka^-.-ak (^i • • • He): 



"-j Ccti (ri ■ ■ - rk) (canonical ensemble) 
J2j exp(/3Aij) Qj g^i-a, (ri, • • • , ffc) (58) 
(grand canonical ensemble). 

With these definitions, and using eqns. (p^, (Q), and (p7|), the free energy 
functional (|^) or (^3|) can be expanded as 

/3.F{p„,c^„} = - ^ jdrSpcir) 5LOp{7^) (59) 

X! ^ X! • ■ ■ ^^fc { (^'l ■ • • rk) Spai (ri) • • • (5pa^ (ffc) 

k = 2 ■ Ql - Qfc 

The field deviations SuJa arc functions of the concentration perturbations by 
means of the relation (|o|) in the canonical ensemble, or the relation ([ij) in the 
grand canonical ensemble. In terms of the correlation functions K introduced 
above, these can be rewritten as 

oo 

5p.{r) = -ETfc^ E K^-'-df, (60) 

irf2 ■ ■ - T^k) SuJa2 {r2) ■ ■ ■ SuJak i^k) ■ 

To leading order in Spa, one gets 

dio^ir) = -Y. fdr" Kj{f, ^Sppif), (61) 

13 

where K~l is the inverse of Kap^ defined through the relation 
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Inserting this into eqn. ( |59| ) yields the gaussian contribution of the fluctua- 
tions to the free energy functional 

f3j: = PF^''^'^ + (62) 

a/3 •' 

The study of these harmonic fluctuations already permits a stability analysis 
of the self-consistent field solution F^'~'^ , and the calculation of scattering 
functions. The partition function in leading harmonic order contains only 
gaussian integrals and can be evaluated in a straightforward way. At first 
sight, inverting the function Kap{f, r') seems a difficult task. However, a con- 
siderable simplification can often be achieved by an appropriate basis trans- 
formation, i.e., expressing Kap, Pa etc. in terms of suitable basis functions 
- for example, in periodically ordered structures an expansion in the Bloch 
waves belonging to the periodically ordered potential uoa (^) has proved useful 
Higher order corrections to eqn. (^2|) can be obtained iteratively, yet 
their treatment becomes much more involved. 

It is often expedient to recast the theory in a way which allows to separate 
total density fiuctuations from composition fiuctuations. For simplicity, we 
shall now consider a mixture which contains only two types of monomers A 
and B of equal size. We replace the variables {pa,Pb) by {p, (f) with the total 
density p = pA+ps, and the composition variable = pa^Pb- The conjugate 
fields are given by tOp — {lua + ujb)/'2 and uj^ — {to a — ws)/2. Furthermore, 
we adopt the definitions of Laradji et al |5J] 



E = Kaa + Kab + Kb A + Kbb 

(63) 



Al = Kaa + Kab — Kb a — Kbb 



A2 = Kaa — Kab + Kba — Kbb 
C = Kaa - Kab ~ Kba + Kbb ■ 

and assume the matrix notation E StOp = jdr'T,{f,r')5L0p{r'). Equation (|6 
can then be rewritten to leading gaussian order 



Sep I Al C j Uco. 



Inverting this expression yields 

(5L0p\ ^ ( E-i A-i 
\5ujJ \ A-i 











= E 


- A2C-IA1 


Al 


= Al 


- CA^^S 


A2 


= A2 


- EA-^C 


C 


= C 


- A1E-1A2 



+ • • • (64) 



(65) 



with X,^^ (66) 
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The free energy expansion (|6^) thus takes the form 

(3T = I3F'^^^ + (67) 
1 S-i+/3V^p^ 




where V^^''(r, r*) is defined as 6'^V/S4'{r) S(j){r') etc., in analogy to eqn. (^. 
Usually in the literature ^ , incompressible melts have been considered 

{dp = 0) with an interaction potential of the form /3Vinter = Jdrx {paPb)/ P- 
In this case, the free energy can be written as 

PJ^ = PF^^P + ]^5p[C''P^]-^5p with C^^^ = [C-i ~ (68) 



Equation (|6^) generalizes this result for arbitrary potentials V and compress- 
ible blends. 

In well segregated macro- or microseparated blends, the main composi- 
tion fluctuations are due to capillary wave fluctuations of interface positions. 
These are less well captured by perturbative treatments, since the local con- 
centration deviations 5pa{r) are not small. As long as the interfaces are still 
localized, i.e., in periodic structures, the method can still be applied to some 
extent [^Sj. However, the expansion will fail if the fluctuations are strong 
enough to destroy the long range order. Note that single, isolated interfaces 
between macroscopic phases usually delocalize ^oj, i.e., the 

extent of the fluctuations depends on the size of the system. In order to study 
these long wavelength fluctuations, other approaches like effective interface de- 
scriptions with input parameters taken from self-consistent field calculations 
are more appropriate. 

In order to illustrate the problem, we shall discuss the simplest example, 
a single interface separating two macroscopic phases in a blend of incompat- 
ible, but otherwise symmetric A and B homopolymers. Let us assume that 
is only slightly distorted from its lowest energy flat state, such that the lo- 
cal deviations of the interfacial position can be described by a single valued 
function h{x,y), and the gradients |V/i| are small. The increase in interfacial 
area caused by these fluctuations cost the free energy 

new = I jdxdy\Vh\\ (69) 

where a is the interfacial tension, and energy losses due to the distortion of 
profiles and the like have been neglected. The functional (|6^) is commonly 
referred to as the capillary wave Hamiltonian. It can be diagonalized by means 
of a Fourier transformation in x and ?/, and since it is quadratic, the spectrum 
and the distribution functions can be determined analytically. The thermal 
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average of the Fourier components is given by = kBT/{(jq^) and one 

obtains a Gaussian height distribution function p5[ 

{5{z - h{x, y))) = (z) = ^ exp(-^) (70) 

with ^' = ^1 dq{\h{q)\') = J-\ni^). (71) 

Here one has to introduce an upper cutoff g^ax and a lower cutoff (j,„i„, since 
the integral J dq/q diverges both at g — > cxd and g — > 0. The lower cutoff is 
obviously given by the system size, (7„-,i„ = 271 /L. An important consequence 
is that the width s of the distribution function P{h) grows logarithmically 
with the system size L, i.e., the interface is marginally rough. 

The value of the upper cutoff is less obvious. Clearly, the capillary wave 
Hamiltonian (^) cannot be expected to provide a good description of poly- 
mer interfaces on all length scales. Ideally, one would hope that one can find 
a microscopic length l/q„,ax beyond which ( |69| ) is valid, and that the sys- 
tem can be studied independently by other means on smaller length scales. 
This implies that the coupling of the long wavelength capillary wave fluc- 
tuations with the local structure on short length scales can be neglected. 
In that case, an approximation makes sense which describes the interfacial 
structure by local "intrinsic" profiles, centered at the local interface position, 
Pa{x,y,z) = pa"*\z — h(x,y)) (convolution approximation [pO) ) . The intrin- 
sic profile characterizes the system on the length scale of 1/ g^^x- When looking 
at the interface on a larger length scale l/qo, one obtains apparent profiles 
p^^^\z), which are broadened by the capillary wave modes with wavevectors 
between qo and q„^^^, 

/OO 1 
dhp^^^'\z-h)P,,,.{h) with = — ln(g_/go). (72) 

Note that since Pg2[z) — J^jlh Pgi2[z — h) P(^g2_gi2^{h), the intrinsic profile 

jOq"*'' (z) can be replaced by pi"^^' (z), and the upper cutoff g^^^ by qo, without 
affecting anything on length scales beyond l/qo. Thus the choice of the upper 
cutoff is largely arbitrary. 

Werner et al^^ have tested this concept by extensive Monte Carlo simu- 
lations of interfaces between immiscible phases in symmetric binary polymer 
blends, within the bond fluctuation model (see Section 3.1). The simulations 
were done in a L x L x D geometry with periodic boundary conditions in 
the L directions, and hard walls in the D direction, which favor one the A 
component and the other the B component. The wall interaction parameters 
were chosen beyond the wetting transition, hence large enough to enforce a 
delocalized AB interface which is on average located in the middle of the 
film. In thin films, the capillary wave fluctuations are limited by the film 
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thickness D rather than by the system size L. However, if the film thickness 
D is chosen large enough compared to L, the interface is essentially free. In 
order to study the interfacial fluctuations, the system was split into columns 
of block size B x B and height D, and the Gibbs dividing surface h{x, y) was 
determined in each column. A typical snapshot of the resulting local interface 
position h{x, y) is shown in Figure ^ for the block size B = 8 lattice constants, 
which is roughly the radius of gyration (Rg « 7) . The monomer profiles were 
then taken relative to this position, and after averaging over all columns, 
the interfacial width w was determined by fitting the order parameter profile 
m{z) = (pAiz) — pb{z)) / p{z)_io a tanh profile, m{z) — mb tanh(z/?i)). Us- 
ing eqn. (|7|), one can show[pi| that the apparent width obtained with this 
procedure is broadened according to 

due to capillary waves, where wq is the intrinsic width and Bq — li^lq^-^^^ 
is the coarse- graining length associated with the upper cutoff q^^^. Figure || 
shows the simulation results for the squared interfacial width as a function 
of block size B in films of lateral size L = 128, for various film thicknesses 
D. They become independent of D for thicknesses larger than D > 48, hence 
the interface can then be considered to be free. For large block sizes B, vo^ 
grows logarithmically with _B, with a slope which is in rough agreement with 
the theoretical prediction 1 / {Aa) (dashed line) . For very small block sizes 
B < A, becomes flat. The self-consistent field value wscf"^ = 20.9 is 
reached at the block size Bq « 7. However, looking just at Figure ^, nothing 
in the shape of the curves indicates that there should be anything special 
about the block size 7 or about = 20.5. The regime of logarithmic growth 
already starts at much smaller block sizes. On the other hand, wscf can be 
made intrinsic width by decree. With the upper cutoff g^^^ defined by the 
corresponding block size Bq, the capillary wave Hamiltonian then provides a 
reasonable description of the simulation data. The question remains which is 
the correct choice of the cutoff Bq. In our example, Bq happens to be exactly 
the gyration radius of the chain, but also twice the statistical segment length 
of a chain {b = 3.05), or one and a half times the interfacial width. Each of 
these is a valid candidate for the cutoff parameter. Semenov[^, for example, 
favors the interfacial width. Here again simulations can help to clarify the 
issue. Recently, Werner et al have repeated the analysis leading to Fig. || 
with systematically varied chain length N and/or the monomer interactions, 
i.e., the x parameter. At fixed = 5.1, they find that the cutoff Bq, 
defined by w(i?o) — wscf, scales like Bq cc VN oc a/T/x cx wssl oc Rg. At 
fixed X — 0.16, Bq first increases strongly with N, but levels off faster than 
Rg oc ViV at the largest chain lengths {N = 256) ||9|, |9|. 

In order to gain a better understanding of this problem, let us go back to 
the case of infinitely long chains and recall that self-consistent field theory is 
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supposed to be exact in this limit. More specifically, it gives the correct free 
energy of the system. On the other hand, the interface has a given, finite sur- 
face tension {e.g., assL as obtained in p2[), hence the interface position will 
fluctuate. Concentration fluctuations may exist in the long chain limit, they 
just do not affect the thermodynamics of the system. However, they influence 
the local structure, such as local concentration profiles or orientational prop- 
erties at interfaces. Self consistent field calculations can thus not be expected 
to give equally good results on all length scales. When used to study local 
structure properties in polymer mixtures, one has to ask for the length scale 
for which the theories describe the system best. The simulations of Werner 
et al suggest that it probably approaches a constant (Bq ^ (3 — 4) wssl) 
in the strong segregation limit (TV — » cxj), but is subject to strong chain end 
corrections. It is worth noting that the cutoff Bq is always much larger than 
the block size at which w'^{B) starts to grow logarithmically as predicted by 
([73I). Up to the length scale of Bq, the structures obtained with self-consistent 
field theories thus average not only over the bulk composition fluctuations, 
but also over the capillary wave fluctuations of the interface position. 



3 Applications 



Self consistent field theories are finding widespread use in numerous con- 
texts of polymer or macromolecular physics. For example, they have been 
employed to calculate complicated phase diagrams of copolymer blends or 
mixtures of copolymers and homopolymers||6^, and to study details of 
density profiles at the internal interfaces in such materials!^, ^ pTf . 

They are frequently applied in surface physics, e.g., in connection with sur- 
face segregation!^ ^ and wetting phenomena^, or for the theoretical 
study of polymers or copolymers which are adsorbed on surfaces or grafted 
to surfaces , 100 101 1. We shall not attempt to give an account of all these 
activities here. Rather, we refer the reader to the various recent reviews on 
the different topics, e.g., Refs. [40, |l], 102|, and illustrate the use of the 



method with two examples which are taken from the present author's own 
research: A comparison between self-consistent field calculations and Monte 
Carlo simulations of polymer conformations at interfaces, and a self-consistent 
field study of the phase behavior of short amphiphilic molecules at surfaces. 



3.1 Conformations of polymers at interfaces 

We consider an interface between coexisting phases in a symmetric mixture 
of homopolymers A and B, in the strong segregation regime {x^ = 17), i.e., 
at temperatures well below the demixing temperature. It is studied by self- 
consistent field theory within the Gaussian chain model and the wormlike 
chain model. 
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The self-consistent field calculations are compared with simulations of the 
bond fluctuation model. The latter represents polymers by chains of spatially 
extended effective monomers, which occupy each a cube of eight neigboring 
sites on a cubic lattice, and which are connected by bonds of length < -\/T0 
lattice spacings. The two types of monomers A and B interact pairwise with 
interactions €aa = ^bb ~ —(-ab ~ ^kBTe, if they are less than lattice 
units apart. The equation of state and the local pair correlations in the bulk 
are well known. Moreover, the model was shown to behave like a dense melt 
at volume fraction 0.5 or density p = 1/16, i.e., the chains have almost ideal 
Gaussian statistics, with known statistical segment length b ~ 3 lattice units. 
The simulations were performed at chain length iV = 32 and interaction 
strength e = 0.1, which corresponds to x = 0.53. 

Figure ^ shows the results for the interfacial profile (|^ and All 
lengths are given in units of wssl = b/y/6x, which is the interfacial width 
in an incompressible mixture of infinitely long polymers [3^, and a "natural" 
unit of length within the self-consistent field theory. In these units, the self- 
consistent field theory predicts that the interfacial thickness decreases with 
increasing chain stiffness r]. However, the statistical segment length b = ay/2rj 
increases in turn, such that the net effect is positive: In absolute units {e.g., 
units of the monomer size a), the interfacial width increases with the chain 
stiffness rj in this regime of relatively small chain stiffness rj. (Note that Morse 
and Fredrickson predict the opposite effect, a decrease of the interfacial width, 
in the limit of large rj. 

We turn to the comparison with the Monte Carlo data. The normalization 
in units of wssl still makes sense, since b and wssl are known, whereas the 
parameter r/ is not. The bare simulation profile is broader by a factor of 
almost two than the self-consistent field prediction (open circles). However, 
this can basically be traced back to the effect of capillary waves, as discussed 
in section 2.4. The simulation profiles obtained after splitting the system into 
blocks of size i? = 8 « i?g as in Ref. [0 (where the much weaker segregated 
case X = 0.16 was studied) are in good agreement with the self-consistent 
field prediction (closed circles) [|62| . Alternatively, one can also correct the 
self-consistent field profile for capillary wave broadening by means of eqn. 

Next, the conformations of polymers in the vicinity of such an interface 
are analyzed. We will examine the conformations of the constituting ho- 
mopolymers, and those of single symmetric A : B diblock copolymers (of same 
length), which adsorb to the interface. The results can again be compared 
with Monte Carlo simulations of very diluted copolymers at a homopolymer 
interface!^. We will limit ourselves to the discussion of the orientational 
properties of the molecules here. 

It is instructive to study separately the orientations of single bonds, of 
chain segments, and of whole chains. In self-consistent field theory, bond 
orientations are conveniently calculated within the wormlike chain model. In 
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the following, we will set 77 — 1/2. At that stiffness, the statistical segment 
length b is identical to the "monomer length" a, hence adjacent monomers 
are essentially uncorrelated. This assumption seems reasonable for the bond 
fluctuation model, as long as no bond potentials have been introduced. (An 
improved guess for 77 would probably be slightly larger than 1/2, since chains 
cannot fold back onto themselves). Note that chains with such a small stiffness 
behave almost like Gaussian chains (cf., e.g., Fig. |5|). 

In order to study profiles of bond orientations 6, we define the bond ori- 
entation parameter 

qir) = — =; , (74) 

which is negative for orientation parallel to the interface, and positive for 
perpendicular orientation. The interfacial profiles of q are shown in Figure 
^. Orientation effects on the bond level are found to be overall very weak. 
Homopolymer bonds tend to orient parallel to the interface in the interfacial 
region. The same holds for most parts of the copolymer; Only very few bonds 
in the central region connecting the A and B block turn perpendicular to the 
interface. For both copolymers (dotted line) and homopolymers (not shown), 
the orientation is weakest at the chain ends. These results are basically in 
agreement with the Monte Carlo data (inset) js^]. 

Further away from the interface, the orientation of copolymer monomers 
is driven by a different effect. The chains are pulled towards the interface 
by one copolymer end, and the bonds align perpendicular as a result. At 
distances of several radii of gyration, the copolymers loose contact to the 
interface and the orientation parameter q drops back to zero. The profile 
of q for copolymers thus reflects two different length scales - a tendency of 
parallel alignment in a region of the extent of the interfacial width, which 
gives rise to the central dip in the profiles, and a force towards perpendicular 
alignment over the length scale of the gyration radius, which is the range of 
the interaction between the copolymer and the interface. Note however that 
most copolymer monomers are located close to the center of the the interface, 
and the net bond orientation is parallel. 

Even though interfaces orient single bonds only very weakly, their effect on 
whole chains is much stronger. The orientation of whole chains involves two 
different factors: The orientation of the gyration tensor at constant total gy- 
ration radius or end-to-end radius, and stretching or compression of the chain 
in one direction. The first factor can not be assessed in a self-consistent field 
calculation with Gaussian or nearly Gaussian chains. However, the simulation 
data shown in Fig. ^js^ indicate that the second effect dominates close to an 
interface: The mean-squared components parallel to the interface xy of the 
end-to-end vectors hardly vary throughout the system, for both homopoly- 
mers and copolymers, whereas they strongly depend on the distance from 
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the interface for the z component. Homopolymers are found to be squeezed 
towards the interface, which leads to an effective parallel orientation. Copoly- 
mers show the inverse behavior, they stretch in the direction perpendicular 
to the interface. The effect is strongest for copolymers centered between one 
and two radii of gyration away from the interface, which are pulled towards 
the interface by their one end, and much weaker for copolymers centered right 
at the interface, which do not feel strong orienting forces. The latter can be 
pictured as consisting of two weakly coupled, almost unperturbed homopoly- 
mer blocks A and B. Indeed, the end-to-end vectors of single blocks centered 
at the interface are oriented parallel to the interface, as shown in Fig. ^. 

To conclude, these examples demonstrate that the self-consistent field the- 
ory allows to study the conformational properties of polymers at interfaces 
in great detail, and that a wealth of information can be obtained from such 
calculations. The predictions of the theory were found to be in overall good 
agreement with Monte Carlo simulations. 



3.2 Amphiphiles at surfaces 

Our second example deals with a somewhat more exotic application of a self- 
consistent fiel d the ory, the study of a coarse-grained model for Langmuir 
monolayers [|65, 105 . These arc monolayers of amphiphilic molecules adsorbed 



on a water surface. If the nonpolar chains of the amphiphiles are sufficiently 
long, one observes experimentally two distinct coexistence regions between 



two dimensional fluid phases on increasing the area per molecule |1 06 1: A 
transition from a highly diluted, "gas" like phase (G) into a more condensed, 
"liquid expanded" (LE) phase, and a second region at higher surface coverage, 
where "liquid condensed" domains are present in a "liquid expanded" environ- 
ment. The coexisting high density phases are true fluid phases, as positional 
correlations within them decay exponentially within a few nanometers. In 
contrast, the directions of the bonds connecting nearest neighbor head groups 



appear to be correlated over tens of micrometers! 107 1, which suggests that 
those phases are probably hexatic. They may be untilted (LS) or tilted (de- 
noted L2 here), with different directions of tilt. The transition between the 
liquid expanded and the liquid condensed phase is the monolayer equivalent 
to the "main transition" in bilayers, where the bilayer thickness jumps discon- 
tinuously as a function of temperature |10^ . The latter is probably relevant 
in biological systems, because it is found in lipid membranes at temperatures 
often close to the body temperature (e.5., 41.5*^(7 in DPPC). 

The problem we wish to address is the following: What is the origin of 
the first order transition between these two fluid phases? The onset of bond 
orientational order cannot account for the discontinuity, since the t rans ition 
between a liquid and a hexatic liquid is of Kosterlitz-Thouless type |l09| and 
thus continuous. There is numerous experimental evidence that the flexibility 
of the polar tails plays a crucial role. In particular, the liquid expanded phase 
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disappears if the chains are made stiff by replacing the hydrogene atoms with 



fluorine atoms] If . The system is thus a good candidate for a self-consistent 
field treatment, which takes due account of the conformational degrees of 
freedom of the chains. 

The amphiphiles were modeled as chains containing one head segment, 
which is confined to a planar surface by a harmonic potential, and seven tail 
segments of diameter Aq and length Iq . The conformational weight is given by 
an expression of the type (||) with an additional bending stiffness contribution 
exp[uU{d)], which favors parallel alignment of adjacent segments {6 = 0). 
The adjustable parameter u determines the stiffness of the chains, and the 
actual form of is not of interest here. Chain segments interact via 

repulsive hard core and long range attractive forces. The self-consistent field 
treatment of the interactions essentially follows the lines of Section 2.3, except 
that in those short chains, one has to account explicitly for the extended 
size the monomers by some appropriate coarse-graining over the center-of- 
mass densities of segments. Furthermore, the interactions have an additional 
anisotropic component, as may result, e.g., from local packing effects. It is 
included perturbatively by adding an orientation dependent term 

PVam = ^0^0 [drdudu p{r, u) p{r, u) v {3{uuf - 1). (75) 

J lOTT 

The parameter v is again adjustable and describes the anisotropy per segment 
of the chains. Within this model, one dimensional self-consistent field profiles 
were calculated in the direction z perpendicular to the surface, hence the 
possibility of lateral order (positional order) was not taken into account. Fig. 
^ shows some typical density profiles. Numerous other quantities were also 
evaluated, such as profiles of the nematic tensor, and the in-plane alignment 
of segments (i| = (ux)'^ + (^!/)^> which is only nonzero when the symmetry of 
the xy plane is broken. The free energy was evaluated according to eqn. (||), 
which allows to calculate phase diagrams by means of a Maxwell construction. 

One finds that such a model indeed displays coexistence between two un- 
tilted fluid phases. The phase behavior is driven by the stiffness u of the 



chains [105 , or alternatively by the anisotropic interactions w||6^. On de- 
creasing the interactions or increasing the chain flexibility, the condensed 
phase merges with the expanded phase in a critical point. On making the 
chains stiffer or the interactions higher, the expanded phase becomes unsta- 
ble. Hence the coexistence region ends in a triple point and a critical point, 
like in experimental Langmuir monolayers (Fig. |l^). In addition, one finds 
a phase with uniform tilt in one direction, which is however metastable and 
buried in the coexistence region (not shown). Since both the effect of segment 
interactions and the chain stiffness go down with increasing temperature, the 
u axis or v axis can be interpreted as temperature axis. 

The nature of the phase transition can thus be analyzed. The expanded 
phase is stabilized by the chain entropy of sufficiently flexible chains. The 
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anisotropic interactions between segments, which have a stronger effect in 
systems stiffer chains, are necessary to bring about a distinct condensed phase. 
The phase transition is driven by the interplay of the entropy of the chains 
and their tendency of paraUel aUgnment. 



Next one may ask what happens to the phase diagram (Fig. 10) if the 
stiffness u is changed substantially. From the previous results, we can infer 
that the transition values for v at the critical and triple point are shifted in 
the opposite direction. In the limit of very stiff chains, the two transitions 
merge, and one is left with only one first order transition, from the gas phase 
directly into the untilted condensed phase. This agrees with the experimental 
results of Ref. |11C]. If the chains are made more flexible, on the other hand, 
the already mentioned tilted phase emerges as additional stable phase: The 
gain of conformational entropy at the corresponding surface densities compen- 
sates in part for the loss of surface energy per chain, and surface coverages 
are stabilized which support collective tilt. This is illustrated in Fig. |ll|. The 
liquid expanded phase can then coexist with either an untilted condensed 
phase or a tilted condensed phase. At even higher chain flexibility, the co- 
existing condensed phase is entirely tilted, and turns into an untilted phase 
via a continuous transition upon further compression of the monolayer (not 
shown). Hence tilted phases are stabilized by chain flexibility. Note however 
that other tilting mechanisms are possible {e.g., resulting from a mismatch 
between head size and chain diameter), which are presumably predominant 
in real monolayers. 

In sum, the self-consistent field analysis of this model lays open a rich 
and complex phenomenology. The different phases at low surface coverage 
are largely recovered. This demonstrates again the power of the method even 
for systems with relatively short chains, even though, unlike in polymers, the 
predictions cannot be expected to be quantitative here. 



4 Outlook 

We have reviewed some recent advances in the self-consistent field approach, 
with the aim to give a flavor of the potential and the limitations of this method 
in the study of complex fluids. Many paths of further developments are possi- 
ble. For example, the combination of self-consistent field theories with other 
more microscopic mean field approaches, like the P-RISM theory [Q, might 
open promising routes to tackle new topical problems like polymer crystal- 
lization. In general, the investigation of interrelations between different mean 
field approaches contributes to a deepened understanding of the individual 



methods! Ill, 112|. Another challenging problem for the future is the for- 
mulation of a general self-consistent field theory for semidilute self-avoiding 
chains, which would bridge between the self-avoiding chain statistics on small 
length scales and the random walk chain statistics on larger length scales 
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(Ref. is a first attempt in this direction). Such a method would allow, for 
example, to study interfaces between hydrophobic and hydrophilic polymer 
components in aqueous environment, which arc of high interest in biology. 
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Figure 1: 

Surface segregation profile of the minority component at coexistence below 
the demixing transition (xc) in a symmetric polymer mixture. Dashed line 
shows the volume fraction profile for tr = 0, where the range of monomer 
interactions is assumed to be zero. The inset shows the total density profile. 
fromRef. Q 
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Figure 2: 

Difference A<i>^ between the volume fraction of the minority component at 
the surface and in the bulk, at two-phase coexistence, vs. l/xNjN = 32). 
Data points show the simulation results of Rouault et al (Ref. (72|): Upper 
bound for A^a (circles), and lower bound (squares) . The arrow indicates the 
value of l/x^ which corresponds to Fig. 
from Ref. |7| 
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Figure 3: 

Typical snapshot picture of the local interface position h{x, y) in Monte Carlo 
simulations of a symmetric homopolymer interface at x = 0.16, N — 32. The 
coarse-graining length i? = 8 is roughly the chain's gyration radius. System 
dimensions are D = 64 and L = 64. 
fromRef. 
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Figure 4: 

Squared interfacial width as a function of block size B in Monte Carlo 
simulations of a homopolymer interface confined in films of various thicknesses 
D with lateral dimension L = 128, at x = 0.16 and N = 32. The horizontal 
dotted line marks the self consistent field prediction for a free interface vo^ = 
20.5, and the vertical dotted line the value of the gyration radius of a chain 
Rg = 7.05. The dashed line indicates the theoretically predicted slope of 
1/(4ct) ln(i?), with a = 0.0156 taken from the self consistent field theory. 
FromRef. Q 
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O Simulation data: apparent profile 
• Profile from block analysis (B=8) 
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Figure 5: 

Concentration profile pa/p* vs. z/wssl (with wssl — b/y/6x) in self- 
consistent field theory for Gaussian chains (solid line) and wormlike chains 
of different chain stiffness 77 (dashed lines), at x = 0.53 and chain length 
N — 32. The results are compared to simulation data without (open circles) 
and with (closed circles) block splitting. The thin dashed dotted line shows 
the profi le fo r ry = 0.5 with capillary wave correction, assuming = 2.5 (see 
eqn. [7^)il0p. 
From Ref. [|62| and §|. 
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Figure 6: 

Orientational order parameter q vs. z/wssl at a homopolymer interface 
(x = 0.53, N = 32) for homopolymer bonds (thick sohd hne) and for bonds 
in an adsorbed diblock copolymer at different positions within the chain: end 
bonds (dotted), bonds linking the two blocks (dashed), bonds in the middle 
of a block (long dashed). The inset shows corresponding Monte Carlo data. 
From Ref. H and |||. 
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Figure 7: 

Mean square end-to-end vector component {Ri^) for i = x,y,z in units of the 
average bulk value vs. the distance z of the center of the end-to-end 

vector from the interface in units of wssl, for homopolymers (long dashed) 
and single copolymer chains (dashed). Parameters are x = 0.53, iV — 32. 
Points are the corresponding Monte Carlo data. The thick solid line gives 
the self-consistent field prediction with capillary wave correction according to 
eqn. |7^, with — 2.5 taken from Fig. |5[ 
FromRef. H 
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Figure 8: 

Mean square end-to-end vector component {i = x, y, z) of the copolymer 

blocks in their minority phase (A block in B phase and vice versa) and in their 
majority phase {A in A, B in S), in units of the average bulk value b'^N/6, 
plotted vs. the distance z of the center of the end-to-end vector from the 
interface in units of wssl, compared to Monte Carlo data. Parameters are 
X = 0.53, N ~ 32. Thick solid line shows the self-consistent field prediction 
corrected for capillary waves according to eqn. [t^, with — 2.5 taken from 
Fig. I 

FromRef. ||. 
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Figure 9: 

Density profiles p in units of 1/ (AqZo) vs. z in units of for different molecular 
areas A/Aq. Long and short dashed lines show the center of mass densities 
of tail and head segments, respectively. Solid line shows the coarse-grained 
density, which accounts for the finite extension of the segments. Parameters 
are u — 2 and v = 13.7. 
From Ref. [||. 
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Figure 10: 

Phase diagram for Langmuir monolayers for stiffer chains {u = 2) in the plane 
of anisotropy per segment v vs. area per molecule A in units of chain diameter 
^0- Dashed area indicates region of two phase coexistence. The gas phase is 
found at much larger values of A. Almost the same phase diagram is obtained 
when varying the stiffness parameter u at fixed anisotropy v — 13.5. (Critical 
point at Uc — 1.95, triple point at ut = 2.05 in arbitrary units) 
FromRef. [||. 
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Figure 11: 

Phase diagram for Langmuir monolayers for flexible chains {u = 1.5) in the 
plane of anisotropy per segment v vs. area per molecule A in units of chain 
diameter Aq . The L2 region corresponds to a phase with broken symmetry in 
the xy plane, where the chains are collectively tilted in one direction. 
FromRef. 11. 
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